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Part I 
The Mathematical Theory of Elasticity 
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Chapter 1 
Mathematical Preliminaries 


In this chapter the basic definitions of vector and tensor algebra, elements of tensor 
differential and integral calculus, and concept of a convolutional product for two 
time-dependent tensor fields are recalled. These concepts are then used to solve 
particular problems related to the Mathematical Preliminaries. 


1.1 Some Formulas in Tensor Algebra 


A vector will be understood as an element of a vector space V. The inner product of 
two vectors u and v from V will be denoted by u - v. If Cartesian coordinates are 
introduced in such a way that the set of vectors {e;} = {e1, e2, e3} with an origin 0 
stands for an orthonormal basis, and if u is a vector and x is a point of E 3 then 
Cartesian coordinates of u and x are given by 


uj=u-e;, xj =x-e (1.1) 
Apart from the direct (vector or tensor) notation we use indicial notation in which 


subscripts range from 1 to 3 and summation convention over repeated subscripts is 
observed. For example, 


3 
u-v=J uni = HiVi (1:2) 


From the definition of an orthonormal basis {e;} it follows that 


e; -ej = ôi; Gj=1,2,3) (1.3) 
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where 6;; is called the Kronecker symbol defined by 


_fi if isj 
igel if ij a4) 


We introduce the permutation symbol £ijķ, also called the alternating symbol, 
defined by 


1 if G, j, k) is an even permutation of(1, 2, 3) 
€ijk = 4 —1 if (i, j, k) is an odd permutation of(1, 2, 3) (1.5) 
0 otherwise, that is, if two subscripts are repeated 


The permutation symbol will be used for the definition of the vector product u x v 
of two vectors u and v 
(u x v); = EijkU j Vk (1.6) 


We may observe that the following identity holds true 
EmisEjks = 5 mj Sik = Ômkôij (1.7) 


An alternative definition of the permutation symbol, given in terms of the vectors 
ej, 18 


Eijk = €; - (ej X ex) (1.8) 
Using this definition of ¢;;,, a generalized form of Eq. (1.7) is obtained 


Sip Sig Sir 
EijkE€par = | jp 9jq jr (1.9) 
Skp Skq Ôkr 


Letting k = r in this identity we obtain Eq. (1.7). 
The permutation symbol ¢;;; can be also used to calculate a3 x 3 determinant 


ay ag a3 
EijkaibjCk = bı bo b3 (1.10) 
ĉi C2 C3 


A second-order tensor is defined as a linear transformation from V to V, that is, a 
tensor T is a linear mapping that associates with each vector v a vector u by 


u = Tv (1.11) 
The components of T are denoted by T;; 
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Tij = ej - Tej (1.12) 


so the relation (1.11) in index notation takes the form 


The Kronecker symbol 6; represents an identity tensor that in direct notation is 
written as 1. 
A product of two tensors A and B is defined by 


(AB)v = A(Bv) for every vector v (1.14) 
Thus, in Cartesian coordinates 
(AB);; = Aix Bkj (1.15) 
The transpose of T, denoted by TT, is defined as a unique tensor satisfying the 
property 
Tu-vy=u-T'v forevery u and v (1.16) 
From this definition it follows that 
Tx = TÄ (1.17) 
If T = TT, then the tensor T is symmetric. Also, if T = —T" then the tensor T is 
skew or asymmetric. Therefore, T is symmetric if T;; = Tj;, and skew if T;; = —Tjj. 


Every tensor T can be expressed by a sum of a symmetric tensor sym T and skew 
tensor skw T, that is, 


T = sym T + skw T (1.18) 
where i 

sym T= >T +T") (1.19) 
and r 

skw T = 5 (T - T’) (1.20) 
In index notation 

Tij = Taj + Tis) (1.21) 

where i 

Taj = 5 (Ty + Tj) (1.22) 
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Tin = 4 Gy = Tj) (1.23) 


Ifa tensor P is skew, then there exists a vector w, called the axial vector corresponding 
to P, such that 


Pu =@ xu for every vector u (1.24) 
and it follows from (1.24) that 
Pij = —EijkOk (1.25) 


and 
@i = —€ijk Pjk (1.26) 


For any tensor T the trace of T is denoted by tr (T). In index notation 
tr (T) =T; (1.27) 
The determinant of T is denoted by det (T) and it is 
Tu T T3 


det (T) = | hi Too T3 (1.28) 
TDi T32 T33 


1.2 Alternative Definitions of a Vector and of a Tensor 
Using an Orthogonal Tensor 


We say that Q is orthogonal if and only if 
Q'Q=QQr=1 (1.29) 
For an orthonormal basis e; and orthogonal Q, the vectors 
e; = Qe; (1.30) 


form an orthonormal basis. 

Also, for two orthonormal bases e; and e; , there exists a unique orthogonal tensor 
Q such that (1.30) holds true. 

For a vector w with components in e; denoted by w;, and with components in e€; 
denoted by w;, we have 


ij 


w. = Qjiw; ., ., (1.31). 
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Wi = ĶjiWi (1.32) 


Similarly, for a tensor T we get 


Ti; = Qei Qj Ti (1.33) 
Tit = Qki Qij Ti; (1.34) 

where 
Qji =e; -ej (1.35) 


The set of three quantities wj, or nine quantities 7;;, referred to e;, which transform 


to another set w; or Tj referred to e}, according to (1.31)-(1.32), or (1.33)-(1.34), 
is defined as a vector, or a tensor, respectively. 


1.3 Further Definitions 


The tensor product of two vectors a and b denoted by a @ b is defined by 


(a ® b)u = (b -u)a for any vector u (1.36) 
In components 
(a @ b);j = aibj (1.37) 
Clearly 
tr(a@b)=a-b (1.38) 


Similarly, the inner product of two tensors A and B is defined by 
A-B=tr (A’B) = AyBy (1.39) 
The magnitude of A is defined by 
|A| = (A- A)! (1.40) 
Also, for any tensor T the following relation holds 
T = Tie; ej (1.41) 
and the nine tensors e; ® ej are orthonormal in the sense 
(e; Q ej) - (ex Q €) = ôkiô ji (1.42) 
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Equation (1.41) constitutes the decomposition formula for T in terms of the nine 


orthonormal tensors e; ® €j. 


Also note that if A is a symmetric tensor and B is a skew tensor then [see MTE-2e, 
Example 2.1.5] 

A-B=0 (1.43) 

and for an orthonormal basis e; and an orthogonal tensor Q the vectors e; = Qe; 


form an orthonormal basis [see MTE-2e, Example 2.1.6], and, for two orthonormal 
bases e; and e’, there exists a unique orthogonal tensor Q such that e; = Qe;. 


1.4 Eigenproblem for a Second Order Tensor 


We call à an eigenvalue corresponding to an eigenvector u of a tensor T if 
Tu = Au (1.44) 


For a symmetric tensor T there exists an orthonormal basis {n;} defined by three 
eigenvectors of T corresponding to three eigenvalues A; of T such that 


3 
T= > i; n; Q n; (1.45) 
i=l 


Here 
Tn; =A; n; (no sum oni) (1.46) 
The inverse of T, denoted by T~!, is defined by 
TT! =T'!T=1 (1.47) 
The tensor T7! is closely related to that of an orthogonal tensor. A tensor A is said 
to be an orthogonal if A is invertible, that is if AT! exists and AT! = AT. Thus, A 
is an orthogonal tensor if and only if 


ATA = AAT =1 (1.48) 


For any invertible tensors A and B 


(AB)! =B-!A7! (1.49) 
(A) =A (1.50) 
(AD! = (A7) (1.51) 
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Also, a tensor A is invertible if and only if its matrix [A] is invertible with [A]~! = 


[A7']. 


1.5 A Fourth-Order Tensor C 


A fourth-order tensor C is defined as a linear transformation that assigns to a second- 
order tensor U another second-order tensor V 


V=C[U] (1.52) 
or, in components, 
Vij = Cijxt Uni (1.53) 
The components C;;,; are defined in terms of the basis {e;} by 
Cijki = (ei Q ej) - C [(e 8 €)] (1.54) 
Let {e} be another orthonormal basis, such that 
e; = Qe; (1.55) 
and let C i jkl be components of C with respect to {e;}. A fourth-order tensor C may 


be also defined as the set of 81 quantities Cmnpq that transform to C i jkl according to 
the formula l 


Ciji = Oni Onj Q pk Qqı Cmnpq (1.56) 
The transpose C! of C is defined as a unique fourth-order tensor that satisfies the 
relation 
A - C [B] = C! [A] - B for all second-order tensors A and B (1.57) 
In components 
Ciri = Cktij (1.58) 
Also, 
|C| = sup {IC[A] |} (1.59) 
|A|=1 


is defined as the magnitude |C| of C. Clearly, 
IC [A]| < |C| |A| for every A (1.60) 
a a gemek ee ee eee An EEA LARE EROR 


For a scalar function f = f(x), x € R CE?, the gradient of f at x is defined by 


u= Vf (1.61) 
or, in components, 
u; =f;=6- Vf (1.62) 
The del operator V is defined by 
ð 
V=e— (1.63) 
OXk 


For a vector function v = v(x), the gradient of v is defined by 


or, in components, 
Vij = (VV)ij = Vij (1.65) 


The divergence of v, div v, and the curl of v, curl v, are defined, respectively, by 


div v(x) = tr(Vv) = vii (1.66) 
and 
(curl v) x a= (Vv — Vv!) a for every a (1.67) 
or 
(curl v); =£ijkVk, j (1.68) 


The symmetric gradient of v, denoted by Vv, is defined by 
F: 1 T 
Vv = sym (Vv) = rial + Vv) (1.69) 


Similarly, if T is a tensor field, the divergence of T and the curl of T are defined, 
respectively, in components, by 


(div T); = Ti, (1.70) 


and 
(curl T);; =€ipg Tjq.p (1.71) 


The Laplacian of a scalar field f, of a vector field v, and of a tensor field T, are 
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Af =f kk (1.72) 
(A v); = Vi,kk (1.73) 
(AT)ij = Tij,kk (1.74) 


Also, instead of A we often use V2 


A = V? (1.75) 


1.7 Integral Theorems 


The divergence theorem for a tensor field T. 


Let T be a tensor field on R C E?, let n be a unit outer normal vector to OR. Then 


fno= faivTav (1.76) 


dR R 


Stokes’ Theorem 
Let u and T denote a vector and tensor fields, respectively, on R, and let C be a closed 
curve in R. Then 


f wsat = f (curtwy-nda (1.77) 
C S 
f Ts dt = J (curl T)'nda (1.78) 
C S 


where S is a surface contained in R and bounded by C, n is the unit vector normal to 
S, and s is a unit vector tangent to C. 


1.8 Irrotational and Solenoidal Fields 


A vector field u on R is said to be irrotational in R if 
curl u = 0 on R (1.79) 


A vector field u on R is said to be solenoidal in R if 
https://gioumeh.com/product/theory-of-elasticity-and-thermal-solution/ 


u -n da = 0 for every closed regular surface S in R (1.80) 


Theorem on irrotational fields 
Let R be a simply connected region of E? 


(a) If u is a vector field on R and curl u = 0, then there exists a scalar field f such 


that u = Vf. 
(b) If T is a tensor field on R and curl T = 0, then there exists a vector field v such 
that T = Vv. 


Theorem on solenoidal fields 


(a) If u is a vector field on R and f u - n da = 0 for every closed surface S C R, 


S 
then there exists a vector field w such that u = curl w. 
(b) If T is a tensor field on R and f T'n da = 0 for every closed surface S C R, 


5 
then there exists a tensor field W such that T = curl W. 


Helmholtz’s Theorem 
If u is a vector field on R then there exist a scalar field f and a vector field v such that 


u = Vf + curl v (1.81) 


and 
div v= 0 (1.82) 


1.9 Time-Dependent Fields 


Let f and g be scalar fields on R x T where R is a region of E*, and T = [0,00) is 
the time interval. The convolution f * g of f and g is defined by 


t 
[f * g], t) = / fa, t—-T) g(x, tT) dt (1.83) 
0 


We list properties of convolution that are useful in applications. 
Let f, g, and h be scalar fields on R x T, continuous in time. Then 
(a) fegaers 
b) (f *g)xh=fx(gxh)=fxgx*h 
C) fxe +h)=f*g+f*h 
(d) fxg=0 => f=0 or g=0 
(e) L{f xg} = L{f} L{g} where L is the Laplace transform with respect to t, that 
is, for any function h = h(x, t) 
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L{h} ca t)dt (1.84) 


0 


The convolution product of two scalar functions can be applied to mixed fields. Thus, 
if u is a vector field and f is a scalar field then 


t 
[f xu] (x, t) = / FS (x, t—T)u(x, tT) dt (1.85) 
0) 


If A and B are time-dependent tensor fields, then 


t 
[A x B](x, t) = [Aw t—t)- B(x, t) dt (1.86) 
0 


If A is a tensor field and u is a vector field, then 


t 


[A * u] (x, tf) = [ae t—tT)u(x, T) dt (1.87) 
0 
In components 
[f xu]; = [f *ui] (1.88) 
[A * B] = [Ajj * Bij] (1.89) 
[A xu]; = [Ajj x uj] (1.90) 


1.10 Problems and Solutions Related to the Mathematical 
Preliminaries 


Problem 1.1. Use the properties of the alternator £; jg introduced in Sect. 1.1 to show 
that 
(a x b) x c = (a - c)b — (b - c)a (1.91) 


(a x b) x (c x d) = [a - (c x d)]b — [b - (c x d)Ja (1.92) 
where a, b, c, and d are arbitrary vectors. 
Solution. To show (1.91) we write the LHS of (1.91) in components and obtain 


[(a x b) x c]; = Eijk (a x b) cx = Eijk Ejpq 4p bq Ck = E jki Ejpq 4p bq cg (1.93) 
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Now, by using the £ — ô identify [see Eq. (1.7)] we obtain 


Ejki Ejpq = Skp Sig = ôkq Sip (1.94) 
Therefore substituting (1.94) into (1.93) we receive 
[(a x b) x c]; = ak ck bi — bk Ck di (1.95) 


or 
[(a x b) x c]; = (a- c)b; — (b - e)a; (1.96) 


Equation (1.96) is equivalent to (1.91), and this proves (1.91). To show (1.92) we 
replace the vector c in (1.91) by c x d, and arrive at (1.92). 


Problem 1.2. Show that for any vector u and a unit vector n the following decom- 
position formula holds true 


u = ut + u! (1.97) 
where 
ut =(u-n)n and ul! =n x (u xn) (1.98) 
Also, show that 
ut - ull! = ; u-n=u'-n, u!-n=0 (1.99) 


Note. If u = u(x) is a vector field defined on a surface S in E?, n = n(x) is a unit 
outward normal vector field on S, and P is a plane tangent to S at x, then ut and ull 
represent the normal and tangent parts of u, respectively, with respect to P. 


Solution. The relation (1.91) in components reads 


u; = (u>); + (ul); (1.100) 
where 
(Ut); = (ug ng) ni, (Ul); = Eijk Ni Ekpquplq (1.101) 
Since by the £ — ô identify 
Eijk Ekpq = Ekij Ekpq = SipSjq — SiqSpj (1.102) 


therefore 


ul); = (Sip Sig — ôiq pj) nj Nq Up 
= Ng Ng Uj — Np Up Nj (1.103) 
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Now, if we substitute (u+); and (ul!); from (1.101) and (1.103), respectively, into 
the RHS of (1.100), and take into account that ng ng = 1, we arrive at the LHS of 
(1.100) which proves (1.100). 

The relations (1.93) hold true, since by (1.101) and (1.103) 


ut -ul = at); (ull); = (ag ng)ni (ui — ni Up np) = (uk ng)? — (uk ny)? = 0 


(1.104) 
ut -n = [(ug ng)nj]nj = ug ny =u: n (1.105) 
and 
ull! -n = (a!) nj; = (ui — ni Up np)ni =0 (1.106) 
Problem 1.3. Show that an alternative form of Eqs. (1.107) and (1.108) in Problem 
1.2 reads 
u =u} +u! (1.107) 
where 
ut =(n@u)n and u! =(1—n@n)u (1.108) 


In Eq. (1.108) the symbol © represents the tensor product of two vectors, and 1 is a 
unit second-order tensor [see Eq. (1.36)]. 


Solution. The tensor product of vectors a and b is defined as a second order tensor 
P with the components 
Pij = aj bj (1.109) 


or in direct notation 
P=a®@b (1.110) 


Therefore, Eq. (1.108) in components read 
(ut); =njujnj and (ul); = (6); — ni nj)uj (1.111) 


Substituting (1.111) into the RHS of (1.107) written in components, we arrive at the 
LHS of (1.107) written in components. This proves (1.107). 


Problem 1.4. Let T = T(x) be a symmetric tensor field defined on a surface S in 
E?,n= n(x) a unit outward normal vector field on S, and P a plane tangent to S at 
x. Show that 

T =T} +T" (1.112) 
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where 
T+ =2sym(n@ Tn) — (n : Tmn Qn (1.113) 
and 
T! = (A-ngmTA-ngQn) (1.114) 
Also, show that 
T+ -T!=0, Tna=Ttn, T'n=0 (1.115) 


Note. The tensors T+ and T!! represent the normal and tangential parts of T, respec- 
tively, with respect to the plane P. 


Solution. Equations (1.113) and (1.114), respectively, in components, take the form 
(TH; = ni Tix nk + nj Tik nk — Nk Tkp Np Mi Nj (1.116) 
and 
(T); = (ôip =ni np)Tpa (qj —Ng Nj) = Tij —nj Tig ng =ni Tpj np+ni nj Np Tp 
(1.117) 
Since T is a symmetric tensor, therefore, 
ni Tik nk = ni Tpj Np (1.118) 
Writing (1.112) in components, and substituting (1.116) and (1.117) into the RHS 
of (1.112) we arrive at the LHS of (1.112) which proves (1.112). 
To prove (1.115); note that 
T+- T! = (T5 "i; (1.119) 
If we note that 


(Sip — ni np)ni =O and (5g; —Ng nj) nj =0 (1.120) 


then substituting (1.116) and (1.117) into (1.119) and taking into account (1.120) we 
obtain (1.115);. To show (1.115)2 we write the RHS of (1.115)2 in components to 


obtain 
(T1n); = (ni Tix nk + nj Tik nk — Nk Tkp Np ni Nj) Nj 
= Tik nk + nj nj ne Tjk — ni Np Nk Tkp 
= Tik ng. (1.121) 
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Finally, writing the LHS of (1.115)3 in components we obtain 


(Tln); = (T); nj = ip — ni np) Tpq (qj — nq ninj (1.122) 


If Eq. (1.120) is substituted into (1.122) we obtain (1.115)3. This completes 
solution to Problem 1.4. 


Problem 1.5. Show that if S is a plane x3 = 0 with the unit outward normal vector 
n = (0, 0, —1) then the decomposition formula (1.107) in Problem 1.4 reads 


Tui Ti2 Ti 0 0 T3 Ti; T 0 
Ti Ta Th3 | = | 0 0 T3}+] Ta Ta 0 
T31 T32 133 Bi 132 133 0 0 0 


Solution. Substituting nı = 0, n2 = 0, n3 = —1, into Eqs. (1.116) and (1.117), 
respectively, in the solution of Problem 1.4, we obtain 


0 0 T3 
T‘~|/0 0 T3 
T31 T32 T33 


and 
Tii T2 0 
T! = | Ti T2 0 
0 0 0 


Hence 
T =T} +T" 


which proves the decomposition formula of Problem 1.5. 


Problem 1.6. Let T be a second-order tensor with components T;;, and let T Æ 0. 


Show that 
det T = £ijkTi1 Tj2 Tha (1.123) 
Epqr (det T) = £ijkTip Tq Tkr (1.124) 
Tip Tig Tir 
EijkE pqr (det T) = Tip Tiq Tir (1.125) 


Tkp Tkq Tkr 
Solution. To show (1.123) we use the result (1.123) of Eq. (1.10) 
a1 a2 43 


Eijk Gi bj ck = | b1 b2 b3 (1.126) 
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where a, b, and ¢ are arbitrary vectors. By letting 


ai = Ti, b =T2, Gc =T; (1.127) 
in (1.126) we obtain 
Ti1 To T31 
&ijk Ti: Tj2 Tk = | T2 T2 T32 (1.128) 
Tis T23 T33 
Since 
det(T) = det(T") (1.129) 


Equation (1.128) is equivalent to (1.123), and this proves (1.123). 
To show (1.124) we let 


qj = Tip, bi = Tig, ci = Tir (1.130) 
in (1.126), where p, q, and r are fixed numbers from the set {1, 2, 3}, and obtain 
Tip Trp T3p 
Eijk Tip Tiq Tkr = Tiq Thq Tq (1.131) 
Tir Tor T3r 


Next, multiplying (1.123) by £pqr, we get 


Epqr det(T) = Epgr £ijk Tit Tj2 Tk3 (1.132) 
Since by Eq. (1.9) 
Spi Spj Ôpk 
Epgr Eijk = | 9qi ôqj gk (1.133) 
Sri Ôrj Örk 


therefore, substituting (1.133) into (1.132) and multiplying T;1, Tj2, and Tk3, respec- 
tively, by the first, second, and third column of the determinant on the RHS of (1.133), 


we obtain 
Ty Tp2 Tp3 


Epgr det (T) = | Ty To T33 (1.134) 
T; T2 T,3 


Since det(T) = det (TT), the RHS of (1.131) is identical to the RHS of (1.134), and 
this proves (1.124). 
Finally, to show (1.125) we multiply (1.124) by £ijx and obtain 


Eijk Epqr (det T) = £ijk Eabe Tap Toq Ter (1.135) 
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or by Eq. (1.9) 


Sia Sib Sic 
Eijk Epgr (det T) = | dja jb Sjc} Tap Toq Ter (1.136) 
Ska Skb Ske 


Now, multiplying the first, second, and third column of the determinant on the RHS 
of (1.136) by Tap, Toq, and Ter, respectively, we get 


Tip Tig Ti 
Eijk Epgr(det T) = | Tip Tjq Tir (1.137) 
Tkp Tkq Tkr 


This proves (1.125), and a solution to Problem 1.6 is complete. 


Problem 1.7. Let T be a second-order tensor with components 7;; such that det 
T Æ 0, and let Î be the tensor with components 


EE | 
Tij = 5 Siva jrs Tpr Tys (1.138) 

Show that . 7 
TT? = T'T = (det T)1 (1.139) 
T! = (det T)'T" (1.140) 


Note. The matrix [Ê j] is called the cofactor of the matrix [T;;], while (7 is called 
the adjoint of [7;;]. 


Solution. The relation (1.139) in components takes the form 
Tx Ti = T} Tj = (det T) ô;j (1.141) 
Using (1.192) we obtain 
Tik T = Tik Tix = Tejab Ekcd Tac Tha Tik (1.142) 
Since, in view of (1.139) in Problem (1.6) 
Epqr (det T") = tie Tpi Taj Trk (1.143) 


and 
det(T") = det(T) (1.144) 


http stheredorsakian (kbh peame writtenip thg fomMast icity-and-thermal-solution/ 


a | 
Tik Th = 58 jab £iap(det T) (1.145) 


Now, using the £ — ô identity (1.7) 
Emis Ejks = Ômj Sik — mk bij (1.146) 


which is equivalent to 
Ejkb Eimb = Smk Oij — Smj Sik (1.147) 


and letting k = m = a in (1.147) we obtain 
Ejab Eiab = 2 8;; (1.148) 
Thus, because of (1.145) and (1.148) we obtain 
Tix Tj; = (det T) 6; (1.149) 
which proves the second part of (1.141). 
To prove that = = 
Tix Ty = Ti Tej (1.150) 
we note that 
ar Pe 1 
Tj Tk; = Tri Tkj = z Pkab fied Tac Toa Tk; 


1 
5 fied &jea(det T) = (det T)5;; (1.151) 


and this completes the proof of (1.139). 
To show (1.140) we note that 


TT !=T'T=1 (1.152) 


and by virtue of (1.139) . 
TT’ = (det T)1 (1.153) 


Multiplying (1.153) by T~!, taking into account (1.152) as well as the relations 
ABC= (AB) C=A (BC) (1.154) 
where A, B, and C are arbitrary matrices, we obtain 
T! = (det TJT! (1.155) 
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Fig. 1.1 Coordinate axes ND 


Problem 1.8. The x; system is obtained by rotating the x; system about the x3 
axis through an angle O0 < 0 < 7/2, as shown in Fig. 1.1. Let T be a symmetric 
second-order tensor referred to the x; system. Show that 


T = 7)! (1.156) 


Tj, = Tii cos” 0 + Tiz sin 20 + Tr sin? 0 


1 
T= z2 — Ti1) sin 20 + Tz cos 20 (1.157) 
Tj = Ty; sin? 0 — T2 sin 20 + Toy cos” 0 


and 


Tiz = T13 cos 6 + Tz sin 0 
Tsz = —T13 sin 0 + Tz3 cos 0 (1.158) 
T33 = T3 


Also, show that an alternative form of the transformation formulas (1.157) and (1.158) 


reads 
Ti + Ty = Ti + Tx 
Thy — Tl, + 2iTi, = exp (210) (Too — Tii + 2i T12) (1.159) 
Tiz — i T33 = exp (i0) (113 — T23) 
T33 = T33 
where i = /—1. 
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nates (r, 0, x3), we find 


Tr + Too = Tii + T2 

Too — Trr + 2iT,ọ = exp(2i0) (T22 — Ty, + 2iT{2) 

T-3 — 1793 = exp (i0) (713 — i723) (1.160) 
T33 = T33 


[Do not sum over r and 6 in Eq. (1.160)]. 
Hint. Use the formula T’ = QTTQ where Q? is the matrix 


cos sind O 
QT =| —sind cos 0 
0 0 1 


Solution. The relations (1.156)—(1.158) follow from the formula 
T=Q'TQ (1.161) 


where QT is the matrix 
cos@ sind 0 
Q' =| —sin@ cos 0 (1.162) 
0 0 1 


To show that (1.156) and (1.158) are equivalent to (1.159) use the identities 


exp(ikO) = cos k0 +isinké, k=1,2 (1.163) 
cos 26 = cos” 6 — sin? 0 (1.164) 
Finally, using the correspondence 


xí =r, x) =0, x3 =z% (1.165) 


and 


Ti = Ty, Toz = To, Tio = T,6 


(1.166) 
Tis = T;3, T3 = 163, T33 = T33 


we transform (1.159) into (1.160), and this completes a solution to Problem 1.8. 


Problem 1.9. A tensor T is said to be positive definite if u - Tu > 0 for every u Æ 0. 
Show that if T is invertible, then TTT and TTT are positive definite. 


Solution. We are to show that TTT and TT T satisfy the inequalities 
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u- (TT )u>0 foreveryu 40 (1.167) 


and 
u-(T' Thu>0 foreveryu 40 (1.168) 


To prove (1.167) and (1.168) note that (1.167) and (1.168) are equivalent to 
(T'u)-(T'u) > 0 foreveryu 40 (1.169) 


and 
(Tu) - (Tu) > 0 for every u 40 (1.170) 


The equivalency is implied by the identities 
u-(TT') u = (Tw) - (Tw) (1.171) 


and 
u-(T! T) u = (Tu) - (Tu) (1.172) 


Now, since T is invertible, T" is invertible. Hence 
Tu 40 foreveryu 40 (1.173) 


and 
T'u 0 for every u Æ 0 (1.174) 


As a result, the inequalities (1.169), (1.170), (1.173), and (1.174) imply that T TT 
and TT T are positive definite, and this completes solution of Problem 1.9. 


Problem 1.10. Show that eigenvalues and eigenvectors for the matrix 


101 
T=|020 (1.175) 
10 3 
are given by 
Ap =2-V2, Wy =2, B= 24V2 (1.176) 
and 
(1) 1 1 1 
n = 
: 4) Lomi) 2) aa. (1.177) 


=0 
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ni e ae (1.178) 
n® =0, nP =, nP =0 


(3) 1 1 1 
ny =+ 
J21+/2./2—/2 
nS) =0 (1.179) 
gl 1 


In Egs. (1.176)—(1.179) A; is an eigenvalue corresponding to the eigenvector n® (i = 
1, 2,3). 


Solution. By using steps leading to a solution of an eigenproblem for the tensor T 
given by (1.175), we find that an eigenvalue à corresponding to an eigenvector n is 
a solution of the algebraic equation 


det (T— 21) =0 (1.180) 

while a unit eigenvector n corresponding to A satisfies the equations 
(T-ADn=0 (1.181) 
In| = 1 (1.182) 


Also, it is easy to check that 41, A2, and A3 given by (1.176) satisfy (1.180); 
and n, n®, and n®) given by (1.194), (1.195), and (1.196), respectively, sat- 
isfy Eqs. (1.181) and (1.182). Hence A; and n® (i = 1, 2, 3) respectively, given by 
(1.176) and (1.177)-(1.179) are the eigenvalues and eigenvectors for the matrix T 
given by Eq. (1.176). This completes a solution to Problem 1.10. 


Problem 1.11. Let T be the tensor represented by the matrix (1.183) in Problem 
1.10, and let {e*} be the orthonormal basis obtained from Eqs. (1.185)—(1.187) in 
Problem 1.10 in which the upper signs are postulated. Define the tensor QT in terms 
of components by 

QT =e -e; (1.183) 


Show that 
T* =Q'TQ (1.184) 


is a tensor represented by a diagonal matrix. Also, compute the components T}, 735, 
and T;5, and show that 
t T*=trT=6 (1.185) 


Solution. The orthonormal basis {ex} obtained from Eqs. (1.185)—(1.187) in Prob- 
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1 1 1 1 1 
= , 0, 
! EE JI EJI V2 -=| 
ež = [0, 1, 0] (1.186) 


€; = 


la a 1 o1 
E 1+2 V2- /2 bya ral 


Since 
e = [1, 0, 0] ’ € = [0, 1, 0], eB = [0, 0, 1] (1.187) 


Eq. (1.183) implies that 


1 1 1 


T T 
— $ =0 

ni pies Je 
of = 1 1 e 

13 /2 24/2 21 
0% =1, 03, =0, (1.188) 
of = 1 1 1 re 
of = 1 1 

33 V2 PER 


Hence, we obtain 
Qı =Q, Q2=0, Q3 = Q} 
Q21 =0, Qn = 0%, Q23=0, (1.189) 
031 = 0%, Q32 = Q05, Q3 = 03; 
Since, by Eq. (1.184), 
Tj; = Qir Tka Qaj (1.190) 


therefore, substituting Tka, Q}, and Qaj from Eqs. (1.183) in Problem 1.10, (1.188), 
and (1.189), respectively, into the RHS of (1.190), we obtain 


2(3 — 22 
Tř = Ba? Ty, = Ti, =0 
P= =e 7 S00 (1.191) 


2 
x x x 
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Equation (1.191) imply that T* is represented by a diagonal matrix. In addition, it 
follows from (1.191) that (1.194) holds true, and this completes a solution to Problem 
1.11. 


Problem 1.12. Prove the following identities in which ¢ is a scalar field, u is a 
vector field, and S is a tensor field on a region R C E*. You may use the £ — 6 


relation. 
curl Vo = 0 (1.192) 
div curl u = 0 (1.193) 
curl curl u = Vdiv u — V7u (1.194) 
curl Vu = 0 (1.195) 
curl (Vu!) = Vcurl u (1.196) 
If Vu = —Vu" then VVu=0 (1.197) 
div curl S = curl div ST (1.198) 
div (curl S) = 0 (1.199) 
(curl curl S)? = curl (curl ST) (1.200) 
curl (ġ 1) = — [curl CADA (1.201) 
div (STu) = u - div S + S - Vu (1.202) 
tr (curl S) = 0 for every symmetric tensor S (1.203) 


If S is symmetric then 
curl curl S = —V°S + 2Ẹ (div S) — VV (tr S) + 1[V7(trS) — div div S] (1.204) 
If S is symmetric and S = G — 1(tr G) then 
curl curl S = —V?G + 2V (div G) — 1 div div G (1.205) 
If S is skew and øo is its axial vector then 
curl S = 1 (div w) — Vw (1.206) 


Solution. To show (1.192) we recall that for any vector field g = g(x) the curl 
operator is defined by [see Eq. (1.68)] 


(curl ø); = Eijk Pk, j (1.207) 
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(curl Vo); = £ijk bok; (1.208) 


